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The time-asymptotic behavior of undamped, nonlinear oscillators with a random frequency is 
investigated analytically and numerically. We find that averaged quantities of physical interest, such 
as the oscillator's mechanical energy, root-mean-square position and velocity, grow algebraically with 
time. The scaling exponents and associated generalized difi'usion constants are calculated when the 
oscillator's potential energy grows as a power of its position: U(x) ~ a;^" for |a::j 00. Correlated 
noise yields anomalous diffusion exponents equal to half the value found for white noise. 
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I. INTRODUCTION 

Randomness in the external conditions entails the parameters of a dynamical system to fluctuate. The extent of 
these fluctuations is independent of any thermodynamic characteristic of the system in contrast to intrinsic fluctuations 
the amplitude of which is proportional to the equilibrium temperature, in accordance with the fluctuation-dissipation 
theorem |^. Usually, external randomness appears as a multiplicative noise in the dynamical equations. The 
interplay of noise and nonlinearity in a system far from equilibrium results in some unusual phenomena In 
fact, the presence of noise dramatically alters the properties of a nonlinear dynamical system both qualitatively and 
quantitatively (for a recent review see Q). For example, it was shown recently that in a spatially-extended system, a 
multiplicative noise, white in space and time, generates an ordered symmetry-breaking state through a non-equilibrium 
phase transition whereas no such transition exists in the absence of noise [g[ |6| . Noise can also induce spatial patterns 
1^, 1^ or improve the performance of a nonlinear device through stochastic resonance Furthermore, even if some 
important qualitative features of a deterministic system survive to external noise, their quantitative characteristics 



may change: a stable fixed point may become unstable |10 , a bifurcation may be delayed (noise-induced stabilization) 
Pl [1^ and scale-invariant properties which manifest themselves as power-laws may be altered with the appearance 
of non-classical scaling exponents . 

The discovery of Brownian motors which are able to rectify random fluctuations into a directed motion (noise- 
induced transport) has triggered a renewed interest in the study of simple one-dimensional mechanical models of 
particles in a potential with random parameters . It is well known that a linear oscillator submitted to parametric 
noise can be unstable even if damping is taken into account ||l^, |l^. This noise-induced energetic instability has 
been observed in diverse experimental contexts such as electronic oscillators [|l6[ 0, nematic liquid crystals and 
surface waves (Faraday instability) ||l^. In engineering fields this instability plays a crucial role in the study of the 
dynamic response of fiexible structures to random environmental loading such as the wave-induced motion of off-shore 
structures or the vibration of tall buildings in a turbulent wind |po|| . The presence of nonlinear friction tends to limit 
the oscillation amplitude: the pendulum with randomly vibrating suspension axis and undergoing nonlinear friction, 
known as the van der Pol oscillator, has been studied in the small noise limit using perturbative expansions ]T6| , pTf . 

In the present work, we consider the motion of an undamped nonlinear oscillator trapped in a general confining 
potential and submitted to parametric random fluctuations. Because there is no dissipation, the energy of the system 
increases with time and we shall show that the position, the momentum and the energy grow as power-laws of time 
with nontrivial exponents that depend on the behavior of the confining potential at infinity | |2^ . A key feature of 
our method is to use the integrability properties of the associated deterministic nonlinear oscillator in order to derive 
exact stochastic equations in action-angle variables. We then use the averaging technique of classical mechanics [ p3[ , 
together with a reduction procedure p3, Ea, to calculate exactly the anomalous scaling exponents, irrespective of the 
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amplitude of the noise. Some of our results were derived before in the particular case of a cubic nonlinearity using an 
energy-enveloppe equation ||26|| . Our method enables us to derive the numerical prefactors appearing in the scaling 
laws (generalized diffusion constants), and our analytical predictions compare very satisfactorily with the numerical 
results. In the case of noise correlated in time, the anomalous diffusion exponents are modified: they can be obtained 
by dimensional analysis arguments and the values thus found also agree with numerical results. Throughout this 
work, crossover phenomena between different scaling regimes are emphasized. 

This article is organized as follows. In section ||, we recall that the energy of a linear oscillator with multiplicative 
noise grows exponentially with time and that this growth may be characterized by a Lyapunov exponent. In section II], 
we analyze the classical Duffing oscillator in presence of parametric noise. Our technique allows us to study precisely 
the long-time behavior of the system. In section |v| , we consider a particle in an arbitrary confining potential that 
grows as a polynomial at large distances. In section we discuss the case of colored noise where the presence of a 
new timescale (the correlation time) leads to a nontrivial crossover from the white noise regime to another scaling 
regime. Our conclusions are presented in section [V^ . In Appendix ^ the nonlinear oscillator in presence of both 
additive and multiplicative noise is briefly studied: we show that at long times the effect of additive noise is irrelevant. 
Appendix ^ is devoted to numerical methods and in Appendix ^ some useful mathematical relations are recalled. 



II. THE LINEAR OSCILLATOR WITH PARAMETRIC NOISE 



In this section we recall known results for an undamped linear oscillator submitted to parametric noise, a generic 
and widely studied model, in order to understand the role of external multiplicative noise |l^, |l^ . The dynamical 
equation for such a system is 

^:,it) + {u;^+m)^{t)=0, (1) 

where x{t) represents the position of the oscillator at time t and to its frequency. The random noise ^(t) is a Gaussian 
white noise of zero mean- value and of amplitude V: 

m) = 0, 

mm) - v5{t^t'). (2) 

The physical interpretation of Eq. (|^) is that the frequency of the oscillator is not constant in time but fluctuates 
around its mean value lo because of randomness in the external conditions (external noise) . When these fluctuations 
are deterministic and periodic in time, Eq. (0) is a Mathieu equation which has been extensively studied p3| . Here, 
we are interested in the case where these fluctuations are random with no deterministic part. The origin, a; = and 
dx/dt = 0, is an unstable stationary solution of Eq. (||). As shown in Ref. this instability can be studied from 
the dynamical evolution of the Probability Distribution Function P(x,v,t) of x and v (with v — x = dx/dt). This 
P.D.F. obeys the Fokker-Planck equation ||^, |^ associated with Eq. (Q): 

dP dP 2 dP V , ^ . 

where Eq. (^ is understood according to Stratonovich rules. 

This Fokker-Planck equation leads to a closed system of ordinary differential equations that couple the n + 1 
moments of order n, i.e., moments of the type {x'^^^v^), where n and k are positive integers and < fc < n: 

= {n-k) (x''-''-^v''+^) - iv^k {x'^-'^+^v''-^) + - k{k - 1) {x^'-'^+^v''-^) . (4) 

The divergence of the moments with time results from the existence of at least one positive eigenvalue of the linear 
system (^). In particular, the mean value of the mechanical energy E of the system {i.e. the sum of its kinetic and 
potential energies) grows exponentially with time: 

{E)^l{v^) + ^c.^x')^e>^\ (5) 
where the growth rate ^ is the positive real root of the equation: 



(6) 
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FIG. 1: Linear oscillator with multiplicative noise: Eq. (|l|) is integrated numerically for D = 1 with a timestep St = 10 . 
Ensemble averages are computed over 10'' realizations. We plot the average {\nE(t)) and the ratio (In E{t)) /t (inset) vs. time 
t. Fig. (a): w = 1, A{uj = 1) = \imt^ (In E{t)) /t = 0.219(3); Fig. (b): degenerate case lo = 0, A(lj = 0) ~ 0.580(5). Both 
estimates of the Lyapunov exponent are in good agreement with the theoretical prediction [E^, Esl. 



It has also been proved that the quenched average of the energy E grows linearly with time, hence the Lyapunov 
exponent A, defined as 

A- hm i(logi?), (7) 

t — >oo t 

is finite and strictly positive |2^, . The positivity of the Lyapunov exponent implies the instability of all moments 
at long times. Note that the growth rate yU, defined in Eq. (||), is larger than the Lyapunov exponent because of the 
convexity inequality, log(£^) > (logE). 

In Figure |^(a) , we present the numerical solution of Eq. (|l|) averaged over a large number of realizations of the noise, 
where the pulsation is w = 1. The algorithm used to solve this stochastic differential equation with multiplicative 
noise is inspired from ||2g[| and explained in Appendix A numerical estimate of the Lyapunov exponent, given in 
(a), agrees very well with the analytic expression of Refs. E8|. The usual statistical equipartition of the 



Fig. |i|(a 

total energy between kinetic and potential contributions is satisfied: (E) = uP' {x^) — {v^}. 

In Figure |l|(b), we show the same quantities for the degenerate linear oscillator obtained by taking lj equal to 0. 
This degenerate case exhibits the same behavior as the generic case and the Lyapunov exponent can be calculated by 



taking the w — > limit in the formulas of Refs. |27, 2^. We conclude that the instability triggered by the noise is the 
dominant effect and that the presence of the linear restoring force —uP'x is irrelevant. 

Hence, in order to avoid an exponential increase of the energy and the amplitude of the oscillator, it is necessary 
to go beyond the linear approximation and consider the effect of nonlinear restoring forces [pO| . 



III. DUFFING OSCILLATOR WITH MULTIPLICATIVE NOISE 



We now analyze the effect of a nonlinear restoring force in an oscillator submitted to an external multiplicative 
noise. In order to preserve the x — > —x symmetry, the nonlinear term has to be odd in the amplitude x. In this 
section, we study the particular case of a cubic nonlinearity 

^x{t) + {u? + i(t)) x(t) + A x(tf = 0. (8) 

The coefficient of the nonlinear term is set equal to unity by rescaling the variable x{t) to x{t) ^/X. The random 
noise ^ is a Gaussian white noise of amplitude V as defined in Eq. (g). The deterministic nonlinear mechanical 
system correponding to Eq. ^ is known as the Duffing oscillator. We shall prove that the cubic term is relevant and 
prevents the average amplitude from growing exponentially. Instead of an exponential behavior, the average energy 
of the oscillator as well as the variances of its position and velocity exhibit a power-law behavior with time. We shall 
calculate exactly the associated scaling exponents. 
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A. The degenerate cubic oscillator 
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FIG. 2: Cubic oscillator: energy equipartition. Eq. (^) is integrated numerically for T) = 1 with a timestep 5t — 510~*. 
Ensemble averages are computed over 10* realizations. The ratios {E{t)) / {v{t)'^) and lv(tY) / {x(t)*') are approximately constant 
and respectively equal to 3/4 and 1. 

The linear part of the restoring force, —oj'^x, is negligible in comparison to the cubic term when the amplitude of 
the oscillator is large. In order to study the long-time behavior of the oscillator, we therefore simplify Eq. (||) to that 
of a degenerate cubic oscillator: 



dt2 



x{t) £,{t) xit) ^xitf = 0. 



(9) 



We first study the deterministic part of Eq. and shall add the noise term afterwards |2^. In one dimension, the 
deterministic cubic oscillator is integrable because the energy E defined as 



E 



-X H — X , 
2 4 ' 



(10) 



is a conserved quantity. The exact solution of the mechanical system, x + x^ — Q, for a fixed value of i?, is given by 



sn I 



dn((4£;)it;^ 
cn({AE)h-^ 



1/2 



(11) 

(12) 



where sn, cn and dn are the classical Jacobi elliptic functions of modulus l/\/2 p2[. The properties of elliptic 
functions that will be useful for our purpose are recalled in Appendix From Eqs' ([ll|) and (|lj), we define the 
action-angle variables of the cubic oscillator ^ . The action variable / corresponds to the area under a constant 
energy curve in phase space: 



[j^^""^^ ^2E-^dx - 8 VT^du^ E^l^ oc ^3/4. 



(13) 



The angle variable (/>, canonically conjugate to the action /, is equal to ^(i?)^ (but for an unimportant additive 
constant) where LoiE) is the frequency corresponding to the energy E. From Eqs. ( pj] ) and (p^), the phase is identified 
as 



0= {\E)-^t. 



(14) 
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This phase variable is defined modulo the common real period AK of the elliptic functions sn, cn, dn that appear in 
Eqs. ( |ll| ) and (p^). The quarter of the period is given by 



X = if(^) = ^/2^^il=^ 1.854. 



(15) 



In terms of the energy-angle coordinates {E, (p), the original variables {x, x) read as follows: 



x = E'^'sd{cb;l/V2) and i = (2i?)V2^!^M^, (16) 

dn (0;1/V2) 

where we have introduced the function sd = sn/dn. Hence, for the free system (without noise) the second-order 
dynamical equation, x + x^ = 0, is equivalent to the following two first-order equations, the first one representing 
energy conservation: 

E = 0, 

^ = {AE)i . (17) 

The presence of external noise spoils the integrability of the dynamical system (|9|) and causes E to grow with time 
by continuously injecting energy into the system. We now take into account the external noise and rewrite the system 
( p^ ) in {E, (f) coordinates. The stochastic evolution equation for the energy is given by 

E = xxm = r^E^" . (18) 

dn3(</>;l/^/2) ^ ' 



-Bd-^^T^-^h (19) 



Inverting the relation (|l^) between and a;, 

and using Eqs. ( [l^ ) and (|8|), we obtain the stochastic evolution of the phase variable: 

With the help of the auxiliary variable Q. defined as 

Vl = 2V2E^/^ , (21) 



we derive a compact and symmetric form for the two stochastic evolution equations (18) and (pOl): 

(1 = ^^%^e(0, (22) 
dn3(0) ^ ' ^ ' 

■ O sd^ 

where the elliptic modulus l/\/2 common to all the elliptic functions has been omitted. We emphasize that the 
coupled equations ( p2[ ) and ( psf ) are mathematically equivalent to the initial system and have been derived without 
any approximation. Moreover, the nature of the parametric perturbation has played no role in the derivation: the 
function ^(t) can be a deterministic or a stochastic function with arbitrary statistical properties. 



We now perform a precise analysis of the long-time behavior of the cubic oscillator driven by a multiplicative 
Gaussian white noise. From a heuristic point of view, we observe from Eq. ( ]2^ ) that Vl undergoes a diffusion process 
and should scale typically as t^/'^. We also notice from Eq. ( p3| ) that, as grows, the phase cj) varies more and 
more rapidly with time. Hence the phase (/> is a fast variable and it is natural to average the dynamics over its 
rapid variations. This averaging process leads to some remarkable and general identities between different physical 
quantitites. Thus, we obtain the average of from Eq. (p^. 



(24) 
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FIG. 3: Cubic oscillator: numerical data, obtained for the same runs as in Fig. is represented by solid lines, and compared 
to the predicted power-law behaviors given by Eqs. ( |46| ) (dashed lines) . We plot the time-dependence of the following averages: 
(a) {E{t)), (b) {x{t)^) and (c) {v{t)^). Data plotted in the insets yield the following estimates of the diffusion constants 
De = 0.031(1), = 0.125(5), D„ = 0.043(3), in excellent agreement with theoretical estimates Eqs. (0), (^) and (^. 



where the overUne indicates the mean-value over a period AK. The coefficient in Eq. (EJ) can be calculated explicitly: 



cn^ (0) 



(25) 



the proof of this identity (and of some forthcoming similar identities) is presented in Appendix ^ Equation (| 
nothing but a statistical equipartition relation between the mean values of the total energy and kinetic energy: 



(26) 



Such an identity provides a generalization of the usual equipartition theorem valid for quadratic degrees of freedom. 
Another statistical equality follows from Eq. ( ^6| ) and the definition (^0|) of the energy E : 

ii') = (^'>- (27) 

We emphasize that these equalities are 'universal' in the sense that they are independent of the form of the noise we 
consider. The derivation of these equipartition relations relies only on the hypothesis that the distribution of the angle 
(j) is uniform over the interval [0,4iir] when t — > oo. In particular, identities (E6|) and (|2^) are valid for multiplicative 
as well as for additive noise (see Appendix]^). Our numerical results (Fig. g) confirm these relations and that the 
averaging procedure is therefore valid. 

The same averaging procedure allows us to derive a closed equation for the stochastic evolution of the slow variable 
ri. We start by writing the Fokker-Planck equation governing the evolution of the P.D.F. Pt(ri, cj)) associated with the 
system ( p^ , p3| ). Since the noise term appears as a multiplicative factor, one must be cautious about the convention 
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FIG. 4: Cubic oscillator. Fig. (a): skewness factors of the position x{t) (dashed line) and velocity v{t) (solid line). Fig. (b): 
flatness factors of the position x{t) (dashed line) and velocity v{t) (solid line). Numerical data is obtained from the same runs 
as in Fig. ^. 



used to define stochastic calculus. Here as well as in the following, we shall use Stratonovich rules because they are 
obtained naturally when white noise is considered as a limit of colored noise with vanishing correlation time ||l|, ^j. 
The Fokker-Planck equation corresponding to Eqs. (^2|) and ( ^ ) reads: 

n 

dtP = -^d^P 
where we have defined: 

m-'-^^^f^ and ,(0)=sd^(^). (29) 
dn ((/)) 

The Fokker-Planck equation (p8| ) written in the variables (SI, (j)) is exact because we study the case of a Gaussian 
white noise. In order to pursue our calculations, we assume that Pi(ri, 4>) becomes independent of </> when t oo, i.e., 
that the probability measure for cf) is uniform over the interval [0,4if]. We now average the Fokker-Planck equation 
(p8|) over the angular variable. We shall use the fact that the average of the derivative of any function is zero: 



d4...) = 0. (30) 

This implies that: 

1^ 



d^, {g{(l))d^g{(l))) = and /(0) = because /(0) = -d,t,g{(f)) . (31) 



Using these properties, and in particular the last identity in (31 ), we derive the phase-averaged Fokker-Planck equation 



d,P=^l^d^P-2dn^j , (32) 
where Pt(r2) is now a function of SI and t only, and where M is given by 

M=7W)-^ r '""'^T'^^^ d0 ^ 2 ^0 ^ ''-\^ 0.3655... (33) 

kJ, dn«(0) f^j^ i5(r(i)r 
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The third equahty is obtained by setting u = sd{(j)) / \/2; the expression in terms of the Euler T function can be found 
in [pit . The averaged Fokker-Planck equation ( P4 ) corresponds to the following effective Langevin dynamics for the 
variable fi, 

n = ^+E{t), (34) 

where the effective noise satisfies the relation 

{Eit)E{t')) ^MVS{t-t'). (35) 

We conclude from Eq. (^4|) that fl exhibits a diffusive behavior. The dynamics of all moments of fl can be deduced 
from Eq. (32), as for example: 

dt{n^) = 3MV, (36) 
dtin-^) - lQMV{n^). (37) 

From Eq. (|36|), we calculate the effective diffusion constant Dfi of fl: 

Dn=3MVc^im6V, (38) 



where (fi^) ~ Dfit when t 00. We conclude from Eq. (37) that has a non-trivial flatness: 

jn^) _ 5 
(n^y ~ 3 • 



(39) 



This last equality shows that although the variable fl diffuses as t^/^, it is not Gaussian (and hence it is not a Brownian 
variable). 

From this statistical information about f2, we now derive the scaling with time of the average of E, a;^, and 
calculate the numerical prefactors (generalized diffusion constants). From Eqs. (]2l]), (|38|) and (|39|), we deduce that: 



^ ' 64 3 64 64 5 (r(i)) 



Using the equipartition relations (|26|) and (27), we find that: 



3^ 5(r(i))« 



8 ^ 16 5(r(i))« 

The prefactor in the mean value of the energy, Eq . ( |4C| ) , agrees with that of |^6| . We have verified numerically the 
power-laws predicted by Eqs. (^0|), (^Tj) and (^) and determined the prefactors. The results are displayed in Fig. ||. 
The numerical values for the exponents and the prefactors are in perfect agreement with the analytical calculations: 
De = 0.03131?^ = 0.125 2) and Dy = G.0417X'2. 

The averaged distribution fonction Pt{fl) can be calculated because Eq. (^2|) is exactly solvable due to its invariance 
under rescalings t — ^ X^t , fl — > Xfl , A being an arbitrary real number (this invariance is the same as that of the heat 
equation). Equation (|3^ ) is solved by using the self-similar Ansatz: 



and the P.D.F of is found to be 



P^{a)^J-n''- —. (43) 
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This distribution function allows us to calculate the P.D.F. in the {x,x) variables in the time-asymptotic regime. In 
particular, the skewness and the flatness factors of the position and of the velocity can be calculated analytically. 
Since both variables x and x are parity symmetric, their skewness vanishes. The flatness is given by 

^^'^ ' ^^'^ - 2.66, (44) 



{v^y 4 dn«(</.) (^24)2 5 ■ • ' 

These values are in excellent agreement with the numerical computations shown in Fig. ^. We notice that the variables 
X and X are non-Gaussian because their flatness differs from the Gaussian value 3. 
To summarize, the following scalings have been derived: 

E - f, 

X ^ t. (46) 

In particular, we conclude that the amplitude a; of a degenerate cubic oscillator does not grow exponentially with 
time but has a normal diffusive behavior. 



B. The Duffing oscillator 



We now study the general case of a non-zero pulsation w: 

^xit) + [uj^ + m) At) + Atf = 0. (47) 

Here the coefHcient of the nonlinear term is taken to be unity and the random noise is Gaussian and white as defined 
in Eq. (|). 

The results of sections || and III A show two regimes: starting from a small initial condition, the amplitude of 
the oscillator grows exponentially with time until x ^ uj where the linear and nonlinear terms are of the same order 
and then the amplitude grows as the square-root of time according to Eq. (^). Because the deterministic system 
corresponding to Eq. (|4^ ) is integrable with the help of elliptic functions, this crossover from exponential to algebraic 
can be derived in a quantitative manner. 

When uj is non zero, Eqs. (|ll| ) and (|l^ ) become: 



1/4 su{{AE + u^)h;k^ 



AE' 



where the elliptic modulus k varies with the energy and is given by 



(48) 



cu({AE + u^)h]k\ 
(2i?)V2 V y ^ (49) 



2^/AE + uo 



We notice that k tends to the limiting value 1 / a/2 when the energy goes to infinity. Defining the angle variable as: 

= sd-fil^±^.;fcV (51) 



'2E 
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FIG. 5: Duffing osciiiator: Eq. (g) is integrated numericaiiy for lo = 1, \ = 10~® and T) = 1 witli a timestep St = 10~*. 
Ensembie averages are computed over 210^ realizations. Fig. (a): {E{t)) vs. t; the dashed fine corresponds to tfie expected 
power iaw {E{t)) ~ _D_b(A) t^. Fig. (b): tlie numericai vaiue of the equipartition ratio {E{t)) / {v{t)'^) changes from a value close 
to 1 to 3/4. Fig. (c): the numerical estimate of the generalized diffusion constant -Db(A) = limt-^ao{E{t)) /t^ is 3.0(4) 10*; the 
dashed line corresponds to the expected diffusion constant -Db(A) = 3.1310*. Fig. (d): {In E{t)) /t vs. time t; the dashed line 
corresponds to the Lyapounov exponent A = 0.22 expected in the linear regime. 



we rewrite the dynamical equation in energy-angle coordinates. However, while deriving the dynamical equation for 
(j) we must remember that the elliptic modulus k depends on the energy E and is, therefore, a function of time. After 
reintroducing the multiplicative noise term, the stochastic Dufhng oscillator in the energy-angle variables becomes: 



A sn(0; fc cn 0; k) 

(4£; + w4)4 dn'^{(^;k) 



6 = UE + u;^)i - .a^ i^E + ^*)sd^i<t>;k) _ / Esd*i<t>;k) _ E snjcj,; k)cn{ct>; k) j;^ sd\0; k)de \ 

{AE + uj^)i \{AE + Lu^)i {AE + u;'^)idii^{(l3;k) ) 

As compared to Eq. (pO|), two supplementary terms appear in Eq. (p3|). These terms are related to dk/dt and are 
proportional to lu'^. 

Although the dynamical equations are more complicated than those of the purely cubic case, the analysis can be 
performed as above. We shall however simplify our discussion here by taking k equal to its asymptotic value l/\/2. 
This approximation isjustified as soon as the energy is large. We also replace the noise ^(t) by the effective noise 
E!(t) defined in Eq. (pq). This second approximation is only qualitatively correct since it amounts to neglecting a 
deterministic force in the effective Langevin dynamics for E. We thus obtain: 

E ^ m ■ (54) 
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We deduce from Eq. (^) that as long as <C w^, the energy behaves as the exponential of a Brownian motion and, 
therefore, increases exponentially with time. However, when E > w^, the nonlinear term becomes important, Eq. js^ ) 
reduces to Eq. (|8|) , and the energy grows as the square of time. 

We expect that the crossover from exponential to algebraic growth will appear when E ^ uj'^, or x ^ lu. Using 
unsealed variables, the balance between linear and non-linear terms in Eq. (|^) is obtained when x — Xc ^ uj/y/X. 
Fig. 1^ demonstrates that the two regimes are observed numerically when the nonlinear coefficient A is very small 
compared to u'^: we use the numerical values u = 1 and A = 10~^. We notice that in the short-time linear regime, 
the usual equipartition relation for a quadratic potential is verified {{E) ~ (i^)), while the exponential growth of the 
energy is characterized by the Lyapunov exponent A defined in Eq. (0). In the long-time regime, the equipartition 
ratio reaches its nonlinear value 3/4, while the energy growth becomes algebraic, with a generalized diffusion constant 



De{X) in good agreement with Eq. (40), up to the expected scahng factor: I?_b(A) = De/\- 



IV. THE GENERAL NONLINEAR OSCILLATOR WITH PARAMETRIC NOISE 

We now consider the case of a particle trapped in a confining potential U (x) and subject to an external noise. This 
mechanical system generalizes the DufRng oscillator studied in the previous section and its dynamics is given by 

l^x{t)+m<t)^^^, (55) 

where £,{t) is the Gaussian white noise of Eq. (||). For the potential to be confining, we must have U +00 when 
|a;| CO. We restrict our analysis to the case where U is a polynomial; in order to respect the x —x symmetry U 
is even in x. Hence, when \x\ 00, 

2n 

U ^ — with n > 2 , (56) 
2n 

the coefficient of x^" is chosen to be l/2n by a suitable rescaling. As before, we expect the amplitude of the oscillator 
to grow without bounds at large times. Keeping only the relevant terms in Eq. ( p5| ) leads to 

^x{i)+m<t)+<tr'-' (57) 

The deterministic version of Eq. (|57|), x + x'^"^^ = 0, is integrable because of energy conservation. For a given value 
of E the deterministic solution is 

X = E^^^"Sn{^{2nE)'^ty (58) 

X = {2n)'^E^^^S',,(^(2nE)'^tj , (59) 



where the energy reads 



E=-x^ + — a;2" . (60) 
2 2n ^ ^ 



The function Sn is defined as the inverse function of an hyperelliptic integral as follows 

SUX)^Y ^X = V^J^ ^f=^-J^l (61) 



1 - liiii 



From this definition, we find a relation between iS„ and its derivative 5^ 



/n \ 2n 
The action-angle variables for the deterministic system are given by 



5M-^ri-^^#^y. (62) 



(63) 



/ = 4y ^J2E^ —dx = 4{2"+'^n)^ E— J Vl^^u^du oc 



^S,^l^]^-y^ 1^^. (64) 



/ X ^ 














Jo 
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The angle variable (j) is well defined modulo the period AKn of the function 5„ where 

When the noise term is taken into account, the energy is not conserved and the stochastic evolution of the variables 
E and cf) is given by 

E = xxm = {"^n)^ Sr,{ct>)S'M) m , (66) 

(2n)^ [2nE) 2„ 

Introducing the auxiliary variable O, 

n = {2n)'^E'^, (68) 
the Eqs. (|66| ) and (|6^) can be written in the simpler form: 

il = (n-l)5„((/.)5;(0)C(i), (69) 

• _ n 5„(0)2 



(27i)Tr 



■m- (70) 



As in the case of the Duffing oscillator, we observe that cj) is a fast variable and therefore we shall average out its 
rapid variations. Starting from (p9|), we write 

/q du Vl - 



(i^) = (2n)"^ 5^(0)2 (i?) = 2 (E) , (71) 

Jo 



the last equality is derived by writing u = Sn((f>), and using Eqs. ( |6^ ) and (|65|). Moreover, the following identity is 
true (as can be shown by integrating — u^" du by parts): 



y2« 



du 



du u2" du , = -n / du v^l - u^" + n / , . (72) 

' VI - u2" io io VI - 



(i^>-^(i^>- (73) 



Substituting this identity in Eq. ([71|) leads to 



From the definition (60) of the energy, we derive another statistical equality: 



(74) 



We emphasize again that these generalized equipartition relations are valid for any type of noise, colored or white, 
additive or multiplicative. The only hypothesis is that the probability distribution function Pi (51, (j)) becomes uniform 
in cf) over the interval [0,4A'„] when t oo. We observe from the numerical simulations presented in Fig. ^.c that 
this condition is very well satisfied. 

The averaging process with respect to the fast variable generates an effective Langevin dynamics for the slow 
variable n{t). Starting from the Fokker-Planck equation for the total P.D.F., Pt(f2, </>), and averaging over leads to 
the following evolution equation for the averaged probability distribution Pf(f2) 



where Mn is given by 



- ^ ^ lo'v^ " « + i r(^)r(^)' ^ ^ 
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FIG. 6: General nonlinear oscillator: Eq. ( |57| ) is integrated numerically for D — 1, with a timestep St, and averaged over 10* 
realizations for n = 2 (2n - 1 = 3), St = 5 10"*; n = 3 (2n - 1 = 5), 5t = 5 10"*; n = 4 (2n - 1 = 7), 5* = 10"". Fig. (a): 
average energy {E{t)) vs. time t. Fig. (b): the limit \imt^oo{E{t)) /t"" , an = yields the following estimates of the diffusion 

constants D^"': D^^' = 0.031(1), D^^' = 0.097(3), D*;** = 0.130(5), in excellent agreement with the predictions of Egs. ^) 
and (^t\). Fig. (c): the measured equipartition ratio {E{t)) / {v{t)'^) is close to the theoretical value given in Eq. ([73|): 3/4 
for n — 2; 2/3 for n = 3; 5/8 for n = 4. Data for = 2 is the same as that of Figs. and H. 



r(.) being the Euler Gamma function. The effective Langevin dynamics for the variable fl is thus 

n-^^ + ^nit), ill) 

where the effective Gaussian white noise S„ [t) satisfies the relation 

{^n{t)^n{t'))^Mn'D5{t~t'). (78) 



The averaged distribution fonction Pt(ri) can be calculated because Eq. ([75|) is exactly solvable due to its invariance 
under rescalings t , f7 \Vl , A being an arbitrary real number. The P.D.F. of J7 is found to be 

1 r OF- 1 



from which we obtain the P.D.F. of the energy 



2(n-l) 
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The long-time behavior of the ampUtude, velocity and energy of the general nonlinear oscillator can now be derived. 
Using the equipartition identity ( fz^ ) and Eq. (|80|), we obtain 

p f 3n+l \ 

{E) ^ , ; {2VMnt)^ oc , (81) 



(2n) — 
2n 



n+1 



(E) oc t— . (82) 



Using Eq (|7|), we find that 



1 f 1 



KnJo Iq7^^ (2n)n(n-l)2 

where we made the change of variables u = 5„(0), and used the following identity (obtained by integrating by parts): 



n + 3 Jo Vl - u^" 



u du 



^Vr^^du=3^ / ^===. (84) 
Finally, we deduce from Eqs. (|58|), ( |80| ) and (|83 



Mn / 2MnVt 



i(2.)ir(^) V(2-)^ 



The analytical results for the nonlinear oscillators with quintic x'^ {n = 3) and hcptic nonlinearity {n — 4) are 
as follows: 

For n = 3, (E) ^ 0.097^^/2 ^ ^^2^ ^ 0.145 1^/^ ^ ^^2^ ^ 0.290 i^/^ . (86) 

For n = 4, (E) = 0.130 1''/^ , (i^) = 0.208 1'^/^ , (x^) = 0.347 1^^^ . (87) 



The formulas (|81|), (^2|) and (|8^) were verified numerically. The scaling exponents and the prefactors given in Eqs. (|8£ 
and (|8^ ) are in excellent agreement with the numerical values, as shown in Fig. ^. 
In conclusion, we have derived the following scaling relations: 



E - t—. 



i2(„-l) _ 



In particular, it should be remarked that x undergoes an anomalous diffusion with time with exponent l/(2n — 2). 
If we make n — > 1 formally, this exponent diverges to infinity: this is consistent with the exponential growth of the 
linear oscillator (see Sec. ||). 



We end this section by considering the case of a general confining potential energy U neither necessarily polynomial 
in X nor even in x. The only requirement is that U +oo when |a;| — s- oo. We discuss the qualitative behavior of E, 
X and X at large times from elementary scaling considerations. Suppose first that U ^ \x\^ for large values of |a;|, r 
being an arbitrary real number. 

• If r > 2, then balance between kinetic and potential energies leads to E ^ x^ ^ x"^; thus the time evolution of 
the energy is given by E' '--^ xi^ ^ E~^2^, From the scaling relations E ^ E/t and ^ ~ t^^^^, we conclude that 

This qualitative argument can be made rigorous by generalizing the results obtained above to non-integer values 
of r. 

• If r < 2, then the potential hi is negligible with respect to the multiplicative noise term and we are back to the 
case of the degenerate linear oscillator. Therefore E, x and x grow exponentially with time. 

If the potential grows exponentially, i.e. lA ^ e^'^ , (3 being a positive real number, then similar considerations lead 
to E t and x ^ t^^"^ (leaving aside logarithmic corrections). We then conjecture that the amplitude x diffuses in a 
logarithmically slow manner: x ~ (Inf)^/''. 
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V. COLORED GAUSSIAN NOISE 

We now consider the case where the Gaussian noise has a non-zero correlation time, and discuss how the previously 
found scalings are modified. The system we want to study satisfies the dynamical equation: 



:x{t) +x{t)r]{t) +x{t) 



2ji-l 



0, 



(89) 



where ry is a colored Gaussian noise of zero mean- value. The statistical properties of rj are determined by 



V 
27 



-\t-t'\/r 



(90) 



where r is the correlation time of the noise. The noise 77 can be obtained from white noise by solving the Ornstein- 
Uhlenbeck equation: 



(91) 



where ^{t) is the Gaussian white noise defined in Eq. 

(a) 



and t, t' » T. 



(b) 




1000 




1000 



FIG. 7: Cubic oscillator with colored noise. Eqs. ( |89| , pl| ) are integrated numerically for n = 2, 2n — 1 = 3, X' = l,r = l, with a 
timestep St — 10~^. Ensemble averages are computed over 10"^ realizations. Fig. (a): the average energy {E{t)) grows linearly 
with time. The inset gives an estimate of the diffusion constant De — 0.047(3). Fig. (b): the time-asymptotic value of the 
equipartition ratio {E(t)) / {v{t)'^) is equal to 3/4. 

Introducing action-angle variables as in Sec. we obtain the same set of coupled Langevin equations, Eqs. ( |69| , [70| ), 
where ^{t) is replaced by the colored noise Ty(t). As emphasized previously, the generalized equipartition relations are 
independent of the nature of the noise: Eqs. (|7^) and (Q) remain valid when the noise is correlated in time. This is 
confirmed by numerical simulations (see Fig. |^?b for n — 2). 

The scalings found in sections III and |^ are deduced by averaging the Fokkcr-Planck equation. Here, we must 
write the evolution equation for the joint P.D.F. of x, v = x and 77, Pt{x, v, rj): 



dP 
'dt 



'V— + (a;^" ^ - x-q)- 



dx 



dv 



1 drjP 
T drj 



V_cPP_ 
2^2 9^ 



(92) 



We perform a scaling analysis of this equation in the spirit of [p4[. Balancing the diffusion term with the time 



derivative leads to 77 ~ t^^^. Then we compare the terms of probability current ^ and 



-xri)P 



balance between these terms is possible only if v 



and 



A consistent 
77. We thus find the following scaling relations 



E 

X 



(93) 
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Thus, we predict that the scaling exponents for colored noise are half the exponents calculated for white noise (|8q). 
Numerical simulations indeed confirm that the average energy of a cubic oscillator {n ~ 2) with colored multiplicative 
noise grows linearly with time (see Fig. |7[a). 

The period T of a deterministic oscillator (without noise) decreases as the energy increases: T ~ iJ^T^r from 
Eq. (|6^). When the equations are written in the energy-angle coordinates, two time scales T and r appear. In the 
regime where t ^ T, the correlation time of the noise is much smaller than the typical variation time of the angle. 
Hence the noise can be considered to be white, the averaging procedure can be applied as in Section |^ and the 
scalings found in (^sh are correct. When T ~ r the noise becomes correlated over a period of the free system and 
can not be treated as white anymore. Now, T ^ t corresponds to i? ~ r "-i which leads to crossover time of the 
order tc^^ ^ r^i^ j.e. ^ t^^. For times larger than t^, the scalings (|93| ) are observed. 



VI. CONCLUSION 



A particle trapped in a confining potential with white multiplicative noise undergoes anomalous diffusion: if the 
confining potential grows as a;^" at infinity, the particle diffuses as D^t^^ . We have calculated the anomalous diffusion 
exponent /?„ = l/(2n — 2), and the coefficient D^- We have found similar laws for the diffusion of velocity and energy. 
Thanks to generalized equipartition identities, we have derived universal relations between the exponents and between 
the prefactors. Our calculations are based on the assumption that in the long-time limit the probability distribution 
function becomes uniform in the phase variable. By averaging out the phase variations, an effective, projected, 
dynamics for the action (or energy) can be defined. This technique enabled us to derive the asymptotic distribution 
law of the energy in the t oo limit, and to calculate its non-Gaussian features (skewness and flatness). Our analytical 
results agree with the numerical simulations within the numerical error bars. Thus, the averaging procedure produces 
very accurate results; it would be an interesting and challenging problem to characterize deviations from our results 
and to calculate subleading corrections. 

In the case of colored multiplicative noise, we have deduced the anomalous diffusion exponents from an elementary 
scaling argument. Our result, supported by numerics, shows that the exponents are halved in presence of time 
correlations. The efficiency of parametric amplification decreases if the noise is coherent over a period of the system 
and therefore the particle diffuses at a much slower rate. In this case, however, the averaging technique is harder 
to apply because the noise itself is averaged out to the leading order. A precise calculation in the case of colored 
Gaussian noise still remains to be done. 

We have considered only Hamiltonian systems, i.e. systems where no friction is present. Nevertheless, if the 
damping is small, the results we have derived for the undamped oscillator remain valid until the crossover time 
(identical to the typical decay time of the energy) is reached. Thegeneral case of a non-linear oscillator with (linear) 
friction leads to interesting results and is currently under study ||35| . 
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APPENDIX A: NONLINEAR OSCILLATOR WITH EXTERNAL AND INTERNAL NOISE 

In this Appendix we discuss the behavior of an oscillator submitted to both additive and multiplicative noises. 
Because we are considering non-dissipative systems, there is no stationary probability distribution; the position, 
velocity and energy of the system satisfy scaling laws. 

We first consider the case where the noise is only additive. The dynamics of a linear oscillator with additive noise 
is exactly solvable and the solution of the equation: 

j-^x{t)+u;^x{t)=m, (Al) 

is given by: 

1 /■* 

x{t) = - sin {uj{t - s)) ^(s) ds . (A2) 
^ Jo 
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(a) 



(b) 
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FIG. 8: Cubic oscillator with additive noise. Eq. (A4) is integrated numerically for V — 1 with a timestep St = 510~*. 
Ensemble averages are computed over 10* realizations. Fig. (a): the average energy {E{t)) grows linearly with time. The inset 
gives an estimate of the diffusion constant when noise is additive: De — 0.505(5). Fig. (b): the time-asymptotic value of the 
equipartition {E{t)) / {v{t)^) is equal to 3/4. 



We conclude from this expression that x{t) grows diffusively with time: 



V ( sin 2ijjt \ 



2uj 



2uj 



(A3) 



Such a diffusive behavior is entirely different from the exponential growth of a linear oscillator subject to parametric 
noise. 

When the nonlinear restoring force is cubic, the linear term can be neglected in the long-time limit and the dynamics 
is given by 



dt^ 



x{t) + x{tf ^ at) ■ 



(A4) 



where £,{t) is the Gaussian white noise defined in Eq. This equation can be analyzed as in Sec. [II A . After 
introducing the energy and angle variables, defined in Eqs. (^|) and (p^), we obtain 



E = xm-m'^ cn(c^;l/x/2) 



2\/2£;i/2 



We can show that the auxiliary variable 8 defined as 



e = E2. 



(A5) 
(A6) 

(AT) 



undergoes a normal diffusion process that can be precisely studied after averaging out the fast variable (j). The P.D.F. 
of the variable 8, the scaling laws for E, x and x, as well as the generalized diffusion constants can be calculated 
exactly, as before. We indicate here only the asymptotic power-laws: 



E ^ t, 
X ^ . 



X 



(A8) 



The scaling exponents for additive noise are different from those for multiplicative noise. For exemple, we observe 
that a particle submitted to an additive noise in a quartic potential is subdiffusive with an anomalous exponent equal 
to 1/2, whereas in the presence of multiplicative noise it behaves diffusively. 
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In the general case, with a nonhnearity of the type x^" ^, with n > 1, one can prove that the variable ©, defined 
in Eq. (A_7), diffuses as t^^"^ in the long time limit. Thus the following scalings are satisfied when t — > oo: 



E 

X 
X 



t, 

1 

1 



(A9) 



All these exponents are smaller than the exponents found for a multiplicative white noise in Eq. (|88|). 



Finally, when both additive and multiplicative noises are present, the oscillator is governed by the equation: 

d^ 

■^X{t) + X{t) UnM{t) + X{tf''-^ = Cadd(t) , (AlO) 

where ^muit and ^add are independent white noises of amplitude I?muit and I?add, respectively. If we study the energy 
variation due to noise, E ~ a:i;^muit + ^Cadd — -E^^muit + -E^^^Cadd, we observe that the first term has a dominant 
effect. From this simple argument, we conclude that the multiplicative noise is expected to dominate over the additive 
noise and therefore, asymptotically, the scaling laws will be those derived for the multiplicative noise alone. However, 
a crossover between the two scalings ( |8^ and ( [A^ ) should be observed by choosing Vmuit ^ T^add- Comparing 
Eqs. ( p8| ) and ([A9|), we find that the effect of the multiplicative noise starts to dominate after a crossover time of the 

order of {Vmuittc)^- - {Vaddtc)^ «.e., ^ V^^^^' /V^^u- 



APPENDIX B: NUMERICAL ALGORITHM 



The algorithm used to integrate numerically the stochastic ordinary differential equations studied in this article 
is the one-step collocation method advocated in ||29|l . In this Appendix, we recall the general principles underlying 
this method, and give the algorithms we used to integrate Eqs. ( |57| ) and (^)-(|9l|), for white and correlated noise 
respectively. All stochastic equations are understood according to Stratonovich rules. 

Let {a;i}i=i,...,Ar be TV real variables of time t, and £^{t) a stochastic process assumed to be Gaussian and white. We 
wish to solve systems of N coupled Langcvin equations of the form 

x^^M{x,{t)})+g,{{x,{t)})at), (Bl) 
where fi and gi are N (smooth) functions of the x^'s. Let St be the integration timestep. Upon integrating formally 



Eq. (Bl) between and St, we obtain the following set of coupled equations implicit in {xi{t)} 



j-St nSt 

x^iSt) - x^{0) ^ M{x,{s)})ds+ g,{{x,{s)})as)ds. (B2) 
Jo Jo 

For St small enough, the functions fi and gi may be Taylor-expanded in the vicinity of t — 0, e.g., 

M{^As)})^M{xj{Q)})+dkM{xj{0)})6xk{s) + ^ {dkdiM{x,{0)})) Sxk{s) 6xi{s) + . . . , (B3) 



where Sxk{s) — Xk{s) — Xk{0). Replacing the functions fi and gi in Eq. (B2) by the expansion (B3), we obtain a new 
set of coupled integral equations implicit in {Sxi}. Upon solving these equations up to an arbitrary order in St, we 
obtain Sxk{St). In practice, we implemented an algorithm exact to 0(5t^). 



1. White noise 



We integrate the following set of first-order differential equations 



X = V, (B4) 
V = -uo'^x- Xx^'^-^ +x£^, (B5) 
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where the Gaussian white noise ^{t) verifies Eq. (^. The exact evolution equations for 5x{St) = x{St) — x{0) and 
Sv{St) = v{St) - v{0) read 

Sx(St) = / ds {v{0) + Sv(s)) , (B6) 







St 



6v{St) = J ds |-cj2 (a;(0) + fe(s))-A (a;(0) + fe(s))^""V(x(0) + (5x(s)) C(s)}. (B7) 
The auxihary variables Zi and Z2 , defined by 



St 



Ziidt) = / as)ds, (B8) 





j-St 

Z2{St) = / ^i(s)ds, (B9) 
Jo 

are Gaussian random variables, with zero average, and the following correlations: {Zi{dt)'^) — Vdt, {Z2{St)'^) — 
{Zi{6t) Z2{5t)) ^V5t^/2. Up to order {5tf , we find 

x{6t) = x{Q)+v{Q)5t + Z2{5t) + ]^6t^ {-Lu'^x{fd)~Xx{Qf''-^), (BIO) 
v{6t) = w(0) + Zi((5t)+(5t (-cj2 2;(0)-Aa;(0)2"-i) +l5t2^(0)(-w2a;(0)-(2n-l)Ax(0)2"^2). (BU) 

In practice, we use two independent Gaussian random noises Zi and Yi with zero mean and correlations (Zi{St)'^) — 
Vdt, (Yi((5t)^) = VSt. As shown in Ref. |29 , the variable Z2 may be approximated by the expression: 



Z2{St) = 6t [ IziiSt) + Yi{St) I , (B12) 



when the algorithm is exact up to order 6t'^. 



2 ' 2V3 



2. Colored noise 



When the noise rjit) is correlated, we must solve a set of three coupled equations 

X ^ V, (B13) 

V = -Lu^x- Xx^'"'-'^ +xr], (BU) 

f] = --Ti+-i. (B15) 
r r 

The set of exact integral equations becomes 

j-St 

6x{St) = / ds {v{0) + Sv{s)) , (B16) 
Jo 

6v{6t) ^ f ds (2;(0) + Sx{s)) - A (x(0) + Ja;(s))^""^ + (a;(0) + Sx{s)) (77(0) + (5r/(s))} , (B17) 



St 



1 . .... 1 



SrjiSt) = ds |-- {rjiO) + 6rjis)) + - • (B18) 

The algorithm to order 6t^ reads 

x{6t) ^ x{0) +v{0)St+^v{Q)6t^ , (B19) 
v{St) ^ v{0) + St {-uj^x{0)- Xx{0)^"-^ +x{0)rj{0)) +^x{0)Z2i5t) 



+ i St^ (^v{0) {-Lj^ x{0) - (2n - 1) A ^(O)^"-^ + .^(q)) ^ i a;(o) r;(0) 



(B20) 



7?(<5i) - ,y(0) + izi(<5t) + -i7?(0)(5t-^Z2((5t) + i,5t2^. (B21) 
r r 2 T-^ 
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APPENDIX C: SOME USEFUL PROPERTIES OF ELLIPTIC FUNCTIONS 



The Jacobian elliptic functions sn(u; fc), cn(u; fc), and dn(u; k) are meromorphic functions of the complex variable u. 
They depend on a parameter k called the elliptic modulus |31, 32l. These functions, like the trigonometric functions, 
have a real period and, like the hyperbolic functions, have an imaginary period. They are thus doubly periodic in the 
complex plane and their periods are {'iK, 2iK'), {AK, 2K + 2iK') and AiK'), respectively, where 



K{k) 



At 



dy 



^{l-t^){l-kH^) Jo y^(i_fc'2y2)(i^_py2) 



with 



y 



1 + A:2y2 



and fc' = 1 - fc2 



(CI) 



The imaginary period is given by K'{k) = K{k') = K (vT^-fc^ • The Jacobian elliptic fimctions satisfy the funda- 
mental relations: 



and their derivatives are 



sv?{(j);k) +cv?{(j}]k) = 1 
fc2sn2(0;fc) +dn^(0;A:) = 1 



(C2) 



— sn((/); fc) = cn((/); fc) dn((/); fc) , 
d0 

— cn((/);fc) = — sn((/); fc) dn(0; fc) , 

-^dn(0;fc) = — fc^sn(0; fc) cn(0; fc) , 



(C3) 



If we define the function sd((/); fc) — sn((/); fc)/dn((/); fc) , the relations (C2) and (|C3|) lead to two new formulas 

^sd(0;fc) = ^4^ and = [l - k'\A\cj,-k)){l + k\d\c^-k)) . 

d^ ' dn2(0;fc) dn^(0;fc) ^ " 



Moreover, from Eqs. ( |C3| ) and (C4), we obtain 

s-ar^{y\k) 
sd-i(2/;fc) 



dt 



^{l-t^){l~kH^) ' 
du 

y^(l-fc'2u2)(l + fc2^t2) 



(C4) 

(C5) 
(C6) 



Similar relations can be derived for cn~^ and dn~^. The expressions of inverse elliptic functions in terms of qu art ic 
integrals provide a classic way of defining elliptic functions. The change of variables t = uj \J\ + k'^v? in Eq. ( |C5| ) 
leads to the identity 



sd "'^(y; fc) = sn ^ 



:;fc 



(C7) 



We deduce from Eq. (C5) 
d 



-^sn ^(y\k) ^ fc 



edt 



= fc 



= fc 



^{\-f^){\-kH-^) Jo ^(l-fc'2^i2)(l + fc2y2) Jo 



sd2(0;fc)d</), (C8) 



where the last equality was obtained by setting u = sd{(/); fc) with = sn ^{y; fc) . 



We verify from Eqs. ( p2[ ) and ( |C3| ) that a; and x defined in Eqs. ( pi] ) and ([T^) are solutions of the cubic oscillator 
X + = 0. Similarly, we verify that x and x defined in Eqs. (|4^) and ( ^9| ) satisfy the differential equation for the 
Duffing oscillator, x + uj'^x + x^ = 0. We now prove the identities used to derive the equipartition relations and the 
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diffusion constants in section [II A. In the case of a purely cubic oscillator, k = k' = l/\/2, and we obtain from 
Eq. Q 



dy 



sd-\y;4=) = ^ 



V2 



(C9) 



From Eq. (C4) we deduce 



cn2(^) 1 



K 



77 / d0- 



''{(j)) J^VT^^du 2 



dn^(<^) KJ, -dn^c^) J^i-^ 



(CIO) 



where we have set u — sd(0)/\/2 to obtain the second equality; the final equality results from Eq. ([72|). Similarly, we 
obtain: 



cn4(0) 1 



^ . , cn4 



(0) _J^{l-u*)idu_i 



"1 dti 



dn«(</>) i^Jo dn«(0) /^^-^ 
The last equality is obtained by the following integration by parts: 

/ l.(l-w'*)idu= / 6u*{l~u*)idu= f 6{l-u^)idu- [ 
Jo Jo Jo Jo 



6(l-u4)5dw; 



thus 



[\l-u'^)idu^- [\l-u^)idu. 
Jo Jo 



Now, using Eq. (CIO), we obtain the numerical value 4/7 in Eq. (pll[). We also have: 



sd2(0) 



1 

KJo 



K 



f 

d<l)sd^{(b) =2-^ 



1 u^du 



f 5 



du 



where A4 was defined in Eq. (|3^). The last equality is obtained by using the following identities 



1 6 J 

u du 



1 ,,3 



2 y/T~l? 2 



Vl-u*dM, 



(Cll) 



(C12) 



u^du 

Vl-u4 



■.du = 



'■VT 



5 rl 

du = - / v^l - u'^du . 
2 Jo 



We end this Appendix by explaining how the equations for the phase dynamics are obtained. Eq. ( pO| ) is derived 
from Eqs. (jcgj) and (|l|): 



d 



dt 



-sd- 



4E 



^^mi-'-^^ip^m. (C13) 



In order to derive the dynamics of the phase of a Dufhng oscillator (|53|), we must differentiate Eq. (|51]) with respect 
to time. We must bear in mind that the elliptic modulus k defined in Eq. ( pO[ ) is now a function of time. Therefore, 
when differentiating the function sd""'^, we obtain two contributions: (i) the derivative with respect to the argument 
produces a resu lt similar to the one for the cubic case, (ii) the derivative with respect to the elliptic modulus, calculated 
using Eqs. (C7) and (28), generates the two last terms in Eq. (^. 
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